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Understanding the H 2 /HI Ratio in Galaxies 
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ABSTRACT 

We revisit the mass ratio ^^i" y between molecular hydrogen (H2) and atomic hydrogen (HI) 
in different galaxies from a phenomenological and theoretical viewpoint. First, the local H2- 
mass function (MF) is estimated from the local CO-luminosity function (LF) of the FCRAO 
Extragalactic CO-Survey, adopting a variable CCM0-H2 conversion fitted to nearby observa- 
tions. This implies an average H2 -density f2 H , = (6.9+2.7)-l(T 5 /!~ 1 and£2 H ,/£2 H i = 0.26+0.11 

in the local Universe. Second, we investigate the correlations between ^^ xy and global 
galaxy properties in a sample of 245 local galaxies. Based on these correlations we intro- 
duce four phenomenological models for ^^" y , which we apply to estimate H2-masses for 

each Hl-galaxy in the HIPASS catalog. The resulting H2-MFS (one for each model for ^^i Xy ) 
are compared to the reference H2-MF derived from the CO-LF, thus allowing us to determine 
the Bayesian evidence of each model and to identify a clear best model, in which, for spi- 
ral galaxies, R^* y negatively correlates with both galaxy Hubble type and total gas mass. 

Third, we derive a theoretical model for ^^i" y for regular galaxies based on an expression for 
their axially symmetric pressure profile dictating the degree of molecularization. This model 
is quantitatively similar to the best phenomenological one at redshift z = 0, and hence repre- 
sents a consistent generalization while providing a physical explanation for the dependence 
of ^^i Xy on global galaxy properties. Applying the best phenomenological model for ^^i Xy 
to the HIPASS sample, we derive the first integral cold gas-MF (HI+H2+helium) of the local 
Universe. 

Key words: ISM: atoms - ISM: molecules - ISM: clouds - radio lines: galaxies. 



1 INTRODUCTION 

The Interstellar Medium (ISM) plays a vital role in galaxies as their 
primordial baryonic component and as fuel or exhaust of stars. Hy- 
drogen constitutes 74 per cent of the mass of the ISM. When it is 
cold and neutral it coexists in the atomic phase (HI) and molecular 
phase (H 2 ). While the former follows a smooth distribution across 
large g alactic substructures, the latt er is found in dense molecular 
clouds (Drapatz & Zi nneckerl 1 1984JJ acting as the sole creches of 
newborn stars. The dissimilar but interlinked roles of HI and H 2 in 
substructure growth and star formation have caused a growing in- 
terest in simultaneous observations of both phases and cosmologi- 
cal simulations that distinguish between HI and H 2 . 

Extragalactic observations of HI often use its prominent 21- 
cm emission line, and currently compri se several thousand galax- 
ies (HI Parkes All Sky Survey HIPASS, iBarnes et alfeOOlh . and a 
maximum redshift of z = 0.2 jVerheiien et alj|2007l) . By contrast, 
most H 2 -estimates must rely on indirect tracers, such as CO-lines, 
with uncertain conversion factors. Consequently, the phase ratio of 
neutral hydrogen R mo \ = dMa 2 /dMui and its value for entire galax- 
ies ^„| Xy s ^h 2 /^hi remain debated, and estimates of the uni- 
versal density ratio ^™™ sc = flt^/^Hi vary by an order of mag- 



nit ude in the local Universe (e.g. 0.14, 0.42, 1.1 stated respec tively 
bv lBoselli et al.ll200llKeres et alfeOOllFukugita et aljl998h . 

Ultimately, the uncertainties of H 2 -measurements hinder the 
reconstruction of cold gas masses M Bas = (Mm + Ma 2 )/fi, where 
/3 ~ 0.74 is the standard fraction of hydrogen in neutral gas with 
the rest consisting of helium (He) and a minor fraction of heav- 
ier elements. The limitations of comparing Mm to M gas caused by 
the measurement uncertainties of Mh, culminate in severe difficul- 
ties to compare statistically tight cold gas-mass functions (MFs) of 
modern cosmological simulations with precise HI-MFs extracted 
from Hi-surveys, such as HIPASS. Both simulations and surveys 
have reached statistical accuracies far better than any current model 
for R^* y , and hence the comparison of observations with simula- 
tions is mainly limited by the uncertainty of R^ xy ■ 

As an illustrati on, Fig. [TJ displays the observed HI-MF from 
the HIPASS sample fewaan etal. 2005) together with several sim- 
ulated HI-MFs. The latter are based on the cold gas masses of 
the simulated galaxies produced by two different g alaxy forma- 
tion models applied to the Mill ennium Simulation dBower et al.l 
120061 : iDe Lucia & Blaizotl feoO?). We have converted these cold 
gas masse s into Hi-masses using four models for R^ l ™ y from the 
literature dYoung & Knezeklll989l ; keres et alJ|20blT fBoselli et al.1 
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Figure 1. Dots represent the observed HI-MF by Zwaan et al. (2005]); lines 
represent sim ulated HI-MFs derived from the semi-a nalytic models by 
iBower et al l 12006, red lines) and|De Lucia & Blaizot (2007, blue lines). 
The fo ur models of l?^f* y were adopte d or derived fr omlYoung & Knezek 
<1989l solid line s), iKeres et all j2003l dashed lines), iBoselli et all feOOl 
dotted lines), and lSautv et alj (2003, dash-dotted lines). 



120021 ; ISautv et alj|2003h . The figure adopts the Hubble constant 
of the Millennium Simulation, i.e. h = 0.73, where h is defined 
by Hq = lOO/ikms -1 MpcT 1 with Hq being the present-day Hub- 
ble constant. The differential gas density 0hi is defined as 0hi = 
dp m /dlogM H i, where Phi(M H i) is the space density (i.e. number 
per volume) of Hi-sources of mass Mm- In Fig. [TJ different models 
for galaxy formation are distinguished by colour, while the models 
of R gaI ™ y are distinguished by line type. Clearly, any conclusion re- 
garding the two galaxy formation models based on their HI-MFs is 
affected by the choice of the model for R^ xy ■ 

This paper presents a state-of-the-art analysis of the galaxy- 
dependent phase ratio , the H2-MF and the integral cold gas- 
MF (HI+H 2 +He), utilizing various observational constraints. In 
Section[2] the determination of H2-masses via CO-lines is revisited 
and an empirical, galaxy-dependent model for the CO-to-H2 con- 
version factor (X -facto r) is derived from direct measurements of 
a few nearby galaxies ( iBoselli et alj|2002l and references therein). 
In Section [3] this model is appli ed to recover an H 2-MF from the 
CO -luminosity function (LF) by Ker es et al. I J2003h . T he resulting 
H2-MF significantly differs from the one obtained bv lKeres et al.1 
(2003) using a constant X-factor. Section [4] presents an indepen- 
dent derivation of the H 2 -MF from a H l-sample with well char- 
acterized sample completeness (HIPASS. lBarnes et alj |2001). This 
approach is less prone to completeness errors, but it premises an es- 
timate of the H 2 /HI-mass ratio R^* y ■ Therefore, we propose four 
phenomenological models of R^ xy (as functions of other galaxy 
properties) and compute their Bayesian evidence by comparing the 
resulting H 2 -MFs to the reference H2-MF derived from the CO-LF. 
This empirical method is supported by Section[5] where we analyt- 
ically derive a galaxy-dependent model for R^* y on the basis of 
the re lation between R mo i and the pressure of the ISM dLerov et"al] 
120081) . A brief discussion and a derivation of an integral cold gas- 
MF (HI+H 2 +He) are presented in Section|6] Section [7] concludes 
the paper with a summary and outlook. 



2 THE VARIABLE CO-TO-H 2 CONVERSION 

2.1 Background: basic mass measurement of HI and H 2 

HI emits rest-frame 1.42 GHz radiation (A = 0.21 m) originating 
from the hyperfi ne spin-spin re laxation. Especially cold HI (T ~ 
50 - 100 K, see lFerrierdl200lh also appears in absorption against 
background continuum sources or other Hi-regions, but makes up 
a negligible fraction in most galaxies. Within this assumption, HI 
can be considered as optically thin on galactic scales, and hence the 
Hl-line intensity is a proportional mass tracer, 



Mh 



2.36 ■ W 



Jykms-' \Mpcj 



(1) 



where 5 hi is the integrated Hl-line flux density and D\ is the lumi- 
nosity distance to the source. 

Unlike Hi-detections, direct detections of H2 in emi ssion rely 
on w eak lines in the infrared and ultraviolet bands dDalgarnol 
2000) and have so far been limited to the Milky Way and a 
few nearby galaxies (e.g. IValentiin & van derWerf 1999). Occa- 
sionally, H 2 has also been detected at high redshift (z x2-A) 
throug h absorptions lines associated with damped Lyman a sys- 
tems dLedoux et alJliooH ; iNoterdaeme et al.ll2008l) . All other H 2 - 
mass estimates use indirect tracers, mostly rotational emission lines 
of carbon monoxide (CO) - the second most abundant molecule in 
the Universe. The most frequently used CO-emission line stems 
from the relaxation of the 7=1 rotational state of the predomi- 
nant isotopomer 12 C ls O. Radiation from this transition is referred 
to as CO(l-0)-radiation and has a rest-frame frequency of 1 15 GHz 
(A = 2.6 ■ 10" 3 m), detectable with millimeter telescopes. The con- 
version between CO(l-0)-radiation and H 2 -masses is very subtle 
and generally expressed by the X-factor, 



X = 



Nn 2 /cm 2 
/co/OCkms- 1 ) 



10- 



(2) 



where Af H2 is the column density of molecules and / C o is the inte- 
grated CO(l-0)-line intensity per unit surface area defined via the 
surface brightness temperature T v in the Rayleigh- Jeans approxi- 
mation. Explicitly, /co = J T v dV = A j T v dv, where V is the radial 
velocity, v is the frequency, and A = |dV/dv| is the wavelength. This 
definition of the X-factor impl ies a mass-luminosity re lation analo- 
gous to Eq. |Q} (see review bv lYoung & Scovilleil99lh . 
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(3) 



where Sco = / Sco.vdV denotes the integrated CO(l-0)-line flux 
and D\ the luminosity distance. Sco.v is the flux density per unit 
frequency, for example expressed in Jy, and thus Sco has units like 
Jykms- 1 . Note that Sco relates to the physical flux F, defined as 
power per unit surface, via a factor A, i.e. F = JSco.vdv = /L'Sco- 
CO-luminosities are often defined as Leo = 4ttD 2 Sco (giving units 
like Jykms-' (hr 1 Mpc) 2 ), thus relating to actual radiative power 
Pco via Aro = 'L'i-co- In the /1-dependent notation above, Eq. (£3j 
remains valid for other molecular emission lines, as long as the X- 
factor is redefined with the respective intensities in the denominator 
of Eq. 10. 



2.2 Variation of the X-factor among galaxies 

The theoretical and observational determination of the X-factor is 
a highly intricate task with a long history, and it is perhaps one of 
the biggest challenges for future CO-surveys. 
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Theoretically, the difficulty to estimate X arises from the in- 
direct mechanism of CO-emission and from the optical thickness 
of CO(l-0)-radiation. CO resides inside molecular clouds along 
with H 2 and acquires rotational excitations from H 2 -CO collisions, 
which can subsequently decay via photon-emission. This mecha- 
nism implies that the CO(l-0)-luminosity per unit molecular mass 
a priori depends on three aspects: (i) the amount of CO per unit 
H2, i.e. the CO/H 2 -mass ratio; (ii) the thermodynamic state vari- 
ables dictating the level populations of CO; (iii) the geometry of 
the molecular region influencing the degree of self-absorption. 

The reason why the CO-luminosity can be used at all as a 
H 2 -mass tracer is a statistical one. In fact, CO-luminosities are 
normally integrated over kiloparsec or larger scales, such as is 
inevitable given the spatial resolution of most extragalactic CO- 
surveys. Therefore, hundreds or thousands of molecular clouds are 
combined into one measurement, and cloud properties, such as ge- 
ometries and thermodynamic state variables, probably tend towards 
a constant average, as long as most lines-of-sight to individual 
clouds do not pass through other clouds, where they would be af- 
fected by self-absorption. The latter as sumption seems correct fo r 
all but nearly edge-on spiral galaxies <Ferrierdl200il :IWall 2006). 
It is hence likely that the different geometries and thermodynamic 
variables of molecular clouds can be neglected in the variations of 
X and we expect X to depend most significantly on the average 
CO/H2-mass ratio of the considered galaxy or galaxy part. How- 
ever, the determination of the CO/IT -ratio is itself difficult and its 
relation to the overall metallicity of the galaxy is uncertain. 

Observational estimations of X require CO-independent H 2 - 
mass measurements, which are limited to the Milky Way and 
a few nearby galaxies. Typical methods use the virial mass of 
giant molecular clouds assum ed to be completely molecular- 
ized < lYoung&Scoville|ll99ll) . the line ratios of different CO- 
isotopomers JWild et al.|[l992h , mm-radiation from cold dust asso- 
ciated with molecular clouds jGuelin et al.ll 19931) , and diffuse high 
energ y y-radiation caused by interactions of cosmic-rays with the 
ISM feertsch et al.lll993l ; lHunter et aljl997t) . 

Early measurements suggested a fairly constant X in the inner 
2-10 kpc of the Galaxy, leading several authors to the conclu- 
sion that X does not signifi cantly depe nd on cloud properties and 
metallicity (e.g. Young & Scoville 1991). This finding has recently 
been supported by Blitz et alJ d2007l) . who analyzed five galaxies in 
the local group a nd found no clear t rend b etwe en metallic i ty and 
X. The results of lYoung & Scovillel jl99ll) and lBlitz et aT] fc007t) 
rely on the assu mption that mo l ecular clouds are virialized. Using 
the same method lArimoto et alj dl996h detected strong variations of 
X amongst galaxies and galactic substructur es, and they f ound the 
empirical power-law relation X oc (O/Hr'. llsraeil 12000) pointed 
out that molecular clouds cannot be considered as virialized struc- 
tures, an d using far-i nfrared measurements rather than the virial 
theorem. Ilsraell 119971) found an even tighter and steeper relation in 
a sample of 14 nearby galaxies, X oc (O/FfT 2 " 7 . 

In summary, despite rigorous efforts to measure X and its rela- 
tion to metallicity, the empirical findings remain uncertain and de- 
pend on the method used to measure X. Since we cannot overcome 
this issue, we shall use a model for X t hat relies on d i fferen t meth- 
ods to measure X, such as presented bv lBoselli et alj 12002'). Their 
sample consists of 14 nearby galaxies covering an order of magni- 
tude in O/H-metalicity. This sample includes early- and late-type 
spiral galaxies, as well as irregular objects and starbursts. For these 
galaxies X was determined from three different methods: the virial 
method, mm-data, and y-ray data. Their data varies from X = 0.88 
in the center of the face-on Sbc-spiral galaxy M 51 to X « 60 
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0.0 


c + ci ■ log(0/H) 


-2.90 ± 0.20 


-1.02 ±0.05 
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c + ci ■ (M B - 5 log h) 


3.67 ± 0.25 


0.176 ±0.006 


0.29 


3.3 


c + ci ■ log(Lco) 


1.85 ±0.15 


-0.288 ± 0.05 


0.29 
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Table 1. Comparison of different models for the X-factor: cq and C[ are 
the best parameters (Gaussian errors are coupled), rms is the rms-deviation 
of the data from the model, and B is the Bayes factor of each model with 
respect to the constant model (first row). 



in NGC 55, a barred irregular galaxy seen edge-on. The high val- 
ues (X > 10) are often associated with dwarf galaxies and nearly 
edge-on spiral galaxies, thus consistent with the interpretation of 
increased CO(1-0) self-absorption in these objects. Typical values 

for non-edge-on galaxies lie around X « 1 - 5. 

For the particular data set of iBoselli et alj d2002l) , we shall 
check the validity of a constant- X model against variable models 
for X, by comparing t heir Bayesian evidence - a powerful tool for 
model selection (e.g. ISivia & SMlingll2006l) . The underlying idea 
is that the probability p{M\d) of a model M given the data set d is 
proportional to the probability p(d\M) of d given M, provided the 
compared models are a priori equally likely (Bayes theorem). The 
probability p(d\M) is also called the Bayesian evidence and can be 
computed as, 



(4) 



p(d\M)= j p{d\0, M)n(e\M)dB 
Jn 

where 6 denotes the vector of free parameters of model M and CI the 
corresponding parameter space; p(d\9, M) designates the probabil- 
ity of the data given a parameter choice 6 and it typically includes 
measurement uncertainties of the data. The prior knowledge on the 
parameters is encoded in the probability density function n(0\M), 
which satisfies the normalization condition j n n(9\M)d9 = 1. Two 
competing models M\ and M 2 are compared by their odds, com- 
monly referred to as the B ayes factor B = p(d\M\)l p(d\M 2 ). Ac- 
cording to Jeffrey's scale ( Jeffrey^ 1 1 96 lh for the strength of evi- 
dence, | In B\ < 1 is inconclusive, while | In B | = 1 reveals positive 
evidence in favour of model Mi (probability=0.750), | In B | = 2.5 
depicts moderate evidence (probability=0.923), and | In B | =5 ex- 
presses strong evidence (probability=0.993). 

We consider the four models listed in Table Q] a constant 
model, where 8 = (co), and three linear models, where 6 = (c , Ci). 
The data are a sample of 14 nearby galaxies, for which X was 
measured (Tabled ; X-factors and O/H-metallicities are taken from 
IBoselli et alj d2002l) and references therei n, while M B -magn itudes 
were taken from the HyperLeda database jPaturel et alj|2003l) . and 
CO(l-0)-luminosities L C o were derived from the references indi- 
cated in Table|2] 

For practical purposes we limit the parameter space CI to 
to e [-10, 10] and c t e [-2,2] and take the prior probabilities 
as homogeneous within CI, i.e. n{0\M) = l/|fl|. The probability 
p{d\0, M) in Eq. ® is calculated as the product, 



p(d\e,M) = [~[ — 1— exp 

1 1 " V27T 



[log(X data ) - log(X™ ode1 )] 2 



2 o- 2 



(5) 



where i labels the different galaxies listed in Tableland cr denotes 
the measurement uncertainty of log(X). We set a equal the average 
value cr = 0.13, for all 14 galaxies. (In fact adopting the specific 
o"- values listed in Table [2] leads to very similar results, but could 
be potentially dangerous as the small value cr = 0.01 of the Milky 
Way is likely underestimated.) 
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Object 


log(0/H) (a) 


,, (b) 
B 

— 5 log h 


log(Lco) 


log(X) (a) 


SMC 


-3.96 


-16.82 


-2.04 (c) 


1.00 


NGC1569 


-3.81 


-15.94 


-1.60 < d > 


1.18 


M31 


-2.99 


-20.23 


-1.40 (e) 


0.38 ±0.21 


IC10 


-3.69 


-15.13 


-1.09^ 


0.82 ±0.12 


LMC 


-3.63 


-17.63 


-0.68 ® 


0.90 


M81 


-3 


-19.90 


-0.07 (h) 


-0.15 


M33 


-3.22 


-18.61 


0.20 (i) 


0.70 ±0.11 


M82 


-3 


-17.30 


0.67 <d) 


0.00 


NGC4565 




-21.74 


1.12 < h » 


0.00 


NGC6946 


-2.94 


-20.12 


124 W 


0.26 


NGC891 




-19.43 


1.48 <h) 


0.18 


M51 


-2.77 


-19.74 


1.80 <h) 


-0.22 


Milky Way 


-3.1 


-19.63 




0.19 ±0.01 


NGC6822 


-3.84 


-16.07 




0.82 ± 0.20 



Table 2. Observational data used for the derivation of a variable X- 
factor (Section 12.2) . Leo is given in u nits of Jykms -1 (tr 1 Mpc) 2 . 
(a) O/H-metallicities and X-factors from iBoselli et all fc002l) . (b) ab- 
solute, extinction-co rrect ed B-Magni t udes f rom the HyperLe d a data base 
<PatureletalJl2003l) . ( c) iRubio et all <199ll) . (d) lYoung et alj Jl989h. (e) 
iHever et alj feOOOh. ffllLerov et alj fe006l) . (g) lFukui et alj Jl999) . (hl lsagel 
ll993l) . (ll lHever et alJ j2004F> ! 



The evidence integrals were solved numerically using a Monte 
Carlo sampling of the parameter space. The resulting Bayes factors 
(listed in TableQJ reveal moderate to strong Bayesian evidence for 
a varia ble X-factor given the X-factors presented by Bos elli et al.1 
(2002). Among the different variable models for log(X), the best 
one depends linearly on log(0/H) (highest Bayes factor), as ex- 
pected from the natural dependence of the CO/H 2 ratio on the 
O/H ratio. However, log(X) is also well correlated with M B and 
log(Lco), and hereafter we will use those relations because of the 
widespread availability of M B and L co data. In fact, a X-factor de- 
pending on L C o simply translates to a non-linear conversion of CO- 
luminosities into H 2 -masses. If the two linear regressions between 
log(X) and M B and between log(X) and log(L C o) were determined 
independently, they would imply a third linear relation between M B 
and log(Lco)- The latter can, however, be determined more accu- 
rately from larger galaxy samples. The sample presented in Section 
14.11 (245 galaxies) yields 

log(Lco) * -4.5 - 0.52 (M B - 5 log h), (6) 

where L C o is taken in units of Jy km s _1 (h -1 Mpc) 2 . To get the best 
result, we imposed this relation, while simultaneously minimizing 
the square deviations of the two regressions between log(X) and 
respectively M B and log(L C o)- In such a way we find 

log(Z) = 1.97-0.308 log(Lco)±ox, (V) 
log(X) = 3.36 + 0.160 (M B - 5 log h) ± cr x . (8) 

These two relations are shown in Fig.[2](red solid lines). For com- 
parison, the independent regressions, obtained without imposing 
the relation given in Eq. 0, are plotted as dashed lines. These rela- 
tions correspond to t he parameters Cp and c t gi ven in Table[TJ Other 
regressions found by Arim oto et alj dl996h and lBoselli et alj j2002h 
are als o displayed. Their approaches are similar, but Arimoto et al.1 
( 1996) used less galaxies (8 instead of 14). The 14 data points in 
Fig.[2]are scattered around the relations of Eqs. and with the 
same rms-deviation of 0.29 in log(X). Combined with the average 
measurement uncertainty of cr = 0.13, this gives an estimated true 
physical scatter in \og(X) of tr x = (0.29 2 - 0.13 2 )" 2 = 0.26. 

The variable models of X given in Eqs. and will be ap- 
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Figure 2. Points represent observed X-factors as a function of CO(1-0)- 
power Lco and absolute blue magnitude Mb for 14 local galaxies. Red 
solid lines represent linear regressions respecting the mutual relation be- 
tween Lco ™d Mb given in Eq. ^6); dashed lines represent indepen- 
dent linear r egressions; the dotted line represents the linear fit found by 
lAnmoto et"alTil996h: and the dash-dotted line represents the linear fit found 



by 



Boselli etalJb002l) . 



plied in Sections[3]and|4] In order to account for the uncertainties of 
X highlighted in the beginning of this section, we shall also present 
the results for a constant X-factor with random scatter in Section[3] 



3 DERIVING THE H 2 -MF FROM THE CO-LF 

Using the variable model for the ^-factor of Eq. 0, we shall now 
recover t he local H 2 -mass fu nction (H 2 -MF) from the CO-LF pre- 
sented bv lKeres et ail (2003). The latter is based on a far infrared- 
selected sub sample of 200 gala xies from the FCRAO Extragalactic 
CO-Survey jYoung et al.1 19951) . which successfully reproduced the 
60/um-LF, t hus limiting the e rrors caused by the incompleteness of 
the sample. Kere s et akl d2003h themselves derived a H 2 -MF using 
a constant model X = 3, which probably leads to an overestimation 
of the H 2 -abundance, especially in the high mass end, where the 
X-factors tend to be lower according to the data shown in Section 

m 

We applied Eq. with scatt er cr x = 0.26 to th e individual 
data points of the CO-LF given bv lKeres et al. The result- 

ing H 2 -MF - hereafter the reference H 2 -MF - is shown in Fig. |3] 
together with the original H2-MF derived by iKeres et al] ( 2003) 
using the constant factor X = 3 without scatter . To both functions 
we fitted a Schechter function dSchechtedl 19761) of the form 



ln(10) • 



Mr: 
M* 



exp 



M* 



(9) 



by minimizin g the weighted squ are deviations of all but the highest 
H 2 -mass bin. lKeresetai] (2003) argue that this bin may contain a 
CO-luminous subpopulation of starbu rst galaxies, sim ilarly to the 
situation in the far infrared continuum ( tYun et alj200ll) . In any case 
the last bin only marginally contributes to the universal H 2 -density. 
The Schechter function parameters are given in Table [3] as well as 
the reduced^ 2 of the fits, total H 2 -densities p H2 an d £2h, = PH 2 /Pcrit» 
and the average molecular ratio ^" 1 "J elsc = Sia^/^m- Both p H2 an d 
Qyi 2 were evaluated from the fitted Schechter function rather than 
the binned data, and Q H i = (2.6 ± 0.3) /r 1 10~ 4 was adopted from 
the HIPASS analysis bv lZwaan et"ai] d2005t) . 

Our new reference H 2 -MF is compressed in the mass-axis 
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Figure 3. Filled squares r epresent our refer ence H2-MF derived directly 
from the observed CO-LF IKeres et al. 2003) using the variable X-factor of 
Eq. j7} with scatter <tx = -26. Open circles are the original H2-MF ob- 
tained bv lKeres et al] J2003l) using a constant factor X = 3 without scatter. 
The dashed and solid lines represent Schechter function fits to our refer- 
ence H2-MF and the original H2-MF, while the dotted line represents the 
Schechter function corresponding to a constant X-factor X = 3 with scatter 
o"x- 

compared to the original one, and our estimate of p H , (Table[3} is 33 
per cent smaller. The global rL/HI-mass ratio drops to 0.26 ±0.11, 
implying a total cold gas density of fl gas = (4.4 ± 0.8) • 10~ 4 h~ l . 
The composition of cold gas becomes: 59 ± 6 per cent HI, 15 ± 6 
per cent H 2 , 26 per cent He and metals, where the uncertainties of 
HI and H 2 are anti-correlated. 

It is interesting to observe the quality of the Schechter function 
fits: the fit to our reference H 2 -MF is much better than the one to the 
original H 2 -MF fares et al.l2003h . Since the original MF is a sim- 
ple shift of the CO-LF (constant .X-factor), the Schechter function 
fit to our reference H 2 -MF is also much better than the fit to the CO- 
LF. We could demonstrate that this difference is partially caused by 
the scatter o"x = 0.26, applied to the variable X-factor when de- 
riving the reference H 2 -MF from the CO-LF. Scatter averages the 
densities in neighboring mass bins, hence smoothing the reference 
MF. Additionally, there is a fundamental reason for the rather poor 
Schechter function fit of the CO-LF: It is formally impossible to de- 
scribe both the H 2 -MF and the CO-LF with Schechter functions, if 
the two are interlinked via the linear transfor mation of Eq. Q. Yet, 
in analogy to the HI-MF fcwaan et al1l2005h . it is likely that the 
H 2 -MF is well matched by a Schechter function, hence implying 
that the CO-LF deviates from a Schechter function. 

We finally note, that the faint end of the reference H2-MF is 
nearly flat (i.e. a = -1), such that the total H 2 -mass is dominated 
by masses close to the Schechter function break at M* ss 10 9 M Q . In 
particular, t he faint end slope i s flatter than for the HI-MF, where 
a = -1.37 iZw aan et alj|2005l) . but it should be emphasized that 
this does not imply that small cold gas masses have a lower molec- 
ular fraction. In fact, the contrary is suggested by the observations 
shown in the Section[4] 

For completeness, we re-derived t he H 2 -MF from th e CO-LF 
using a constant X-factor X = 3 (like lKeres et al.|[2003h with the 
same Gaussian scatter tr x = 0.26 as used for our variable model 
of X. The best Schechter fit for the resulting H 2 -MF is also dis- 
played in Fig. [3] The difference between this H2-MF and the orig- 





reference H2-MF 
(variable X) 


original tb-MF 
(constant X) 


M* 


7.5- 10 8 /T 2 M Q 


2.81 ■ 10 9 /T 2 M G 


a 


-1.07 


-1.18 


</>' 


0.0243 h 3 Mpc~ 3 


0.0089 h 3 Mpc~ 3 


Red.^ 2 


0.05 


2.55 


PH, 


(1.9 ±0.7)- 10 7 AM o Mpc- 3 


(2.8± 1.1)- 10 7 h M Mpc" 3 




(0.69 ± 0.27) ■ 10~ 4 ft -1 


(1.02 ±0.39)- 10- 4 /i-' 


^universe 
mol 


0.26 ±0.11 


0.39 ±0.16 



Table 3. Schechter function parameters, reduced x 1 , and universal mass 
densities as obtained by integrating the Schechter functions. ^™[' erse = 
Cl^/Ckn is the global Fb/HI-mass ratio of the local Universe. The very 
small reduced^ 2 of our reference H2-MF arises from a spurious smoothing 
introduced by the scatter o"x- 



inal H 2 -MF by IKeres et alj d2003h demonstrates that the scatter of 
X stretches the high mass end towards higher masses. 



4 PHENOMENOLOGICAL MODELS FOR THE 
Hj/HI-MASS RATIO 

In this section, we shall introduce four phenomenological models 
for the H2/HI-mass ratio R^* y of individual galaxies. Each model 
will be used to re cover a H 2 -MF from the HIPASS Hi-catalog 
l lBarnes et alj|200lh . thus demonstrating an alternative way to de- 
termine the H2-MF to the CO-based approach. Comparing the H2- 
MFs of this section with the reference H 2 -MF derived from the CO- 
LF (Section[3) will allow us to determine the statistical evidence of 
the models for i? gJl ? xy . 

mol 

4.1 Observed sample 

The sample of galaxies used in this section is presented in Ap- 
pendix [A] and consists of 245 distinct objects with simultaneous 
measurements of integrated Hl-line fluxes and CO(l-0)-line fluxes. 
The latter were drawn from 9 catalogs in the literature, and, where 
not given explicitly, recomputed from indicated H 2 -masses by fac- 
toring out the different X-factors used by the authors. HI line fluxes 
were t aken from H IPASS via the optical cross-match catalog HOP- 
CAT dPovle et al. 2005). Both line fluxes were homogenized using 
^-dependent units, where they depend on the Hubble parameter h. 
Additional galaxy properties wer e adopted from the homogenous 
reference database "HyperLeda" jPatu rel et al. 2003). These prop- 
erties include numerical Hubble types T, extinction corrected blue 
magnitudes Mb, and comoving distances D\ corrected for Virgo in- 
fall. In the few cases, where these properties were unavailable in 
the reference catalog, they were copied from the original reference 
for CO-fluxes. For each galaxy we calculated HI- and ft-masses 
using respectively Eqs. (TJ and (3). The variable X-factors were de- 
termined from the blue magnitudes according to Eq. {8}. We chose 
to compute X from Mb rather than from Leo, because of the smaller 
measurement uncertainties of the M B data. Finally, total cold gas 
masses M gas = (M ffl + M H ,)/y3 and mass ratios = M H ,/M H i 

were calculated for each object. While the masses depend on the 
distances and hence on the Hubble parameter h, the mass ratios 
R^ y = Mn 2 /M m are independent of h. 

This sample covers a wide range of galaxy Hubble types, 
masses, and environments, and has 49 per cent overlap with the 
subsample of the FCRAO Extragalactic CO-Survey used for the 
derivation of the reference H2-MF in Section [3] We deliberately 
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Figure 4. H2/HI-mass ratio versus numerical Hubble type T. Grey dots 
represent the empirical data obtained by applying the variable X-factor of 
Eq. {8} with scatter to the CO-measurements. Black points represent the 
binned data; vertical bars represent statistical uncertainties obtained via 
bootstrapping, i.e. they depict a l-cr confidence interval of the bin average 
obtained by examining 10 4 random half-subsets of the full data; horizontal 
bars represent the bin intervals. The solid line represents model 1 fitted to 
the data points. Open circles and dashed b ars denote the binned data of the 
original paper by Young & Kne zekl ll989t) . 



limited the sample overlap to 50 per cent in order to control possi- 
ble sample biases. 

We emphasize that this sample exhibits unknown complete- 
ness properties, which a priori presents a problem for any empirical 
model for ^^„" y - However, as long as a proposed model is formally 
complete in the sense that it embodies the essential correlations 
with a set of free parameters, these parameters can be determined 
accurately even with an incomplete set of data points. The difficulty 
in the present case is that no reliable complete model for the molec- 
ular fraction R^™ y has yet been established. We shall bypass this 
issue by proposing several models for R^'"^ that will be verified 



with hindsight (Section |4.2t . Additional verification will become 
possible in Section [5] where we shall derive a physical model for 



R: 



.galaxy 



galaxy 



4.2 Phenomenological models for R\ 

The galaxy sample of Section 14.11 reveals moderate correlations 
between ^^ xy and respectively T, M eas and Mb- These correla- 
tions motivate the models proposed below. Other correlations were 
looked at, such as a correlation between R gd] '" iy and environment, 

' mol ' 

which may be suspected from stripping mechanisms acting differ- 
ently on HI and H2. However no conclusive trends could be iden- 
tified given the observational scatter of 7?j?!f y . All our models are 
first presented with free parameters, which are fitted to the data at 
the end of this section. 

such 

— mm 

as is often used in the literature, 



Model (tf^icf) assumes a constant H 2 /HI-ratio R^ y 



lo g(Cro y ) = 9o+^phy.o, (10) 

where q is a constant and <r phy denotes an estimate of the physical 
scatter of perfectly measured data relative to the model. 

Model 1 is galaxy-type dependent, as suggested by earlier 
studies revealing a trend for R^ xy to increase from late-type 
spiral galaxies t o early-type ones (e.g. lYoung & Knezikl 1 19891 : 
ISautv et al]2003l) . The type dependence of our sample is displayed 
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Figure 5. H2/HI-mass ratio versus total cold gas mass M eas = (Mm + 
Mh 2 )/j3. Grey dots represent the empirical data obtained by applying the 
variable X-factor of Eq. (8) with scatter to the CO-measurements. Black 
points represent the binned data; vertical bars represent the l-cr confidence 
intervals; horizontal bars represent the bin intervals. The solid line repre- 
sents model 2 fitted to the data points. 



in Fig. [4] The binned data clearly show a monotonic increase of the 
molecular fraction by roughly an order of magnitude when pass- 
ing from late-type spiral galaxies (Scd-Sd) to early-type spiral and 
lenticular galaxies (SO-SO/a). The unbinned data illustrate the im- 
portance of parameterizing the physical scatter. The Hubble type 
dependence can be widely explained by the effect of the bulge com- 
ponent on the disc size, as detail ed in Section [5] Observa tionally, 
this dependence was first noted by Young & Knezek ( 1989), whose 
bins are also displayed in the figure. Their molecular fractions are 
generally higher, partly due to their rather high assumed X-factor 
of 2.8. The monotonic trend seems to break down between lentic- 
ular and elliptical galaxies, where the physical situation becomes 
more complex. In fact, many elliptical galaxies have molecular gas 
in their center with no detectable Hi-c ounterpart, while others s eem 
to have almost no H 2 (e.g. M 87, see iBraine & W iklind 199J), or 
may even exhibit Hi- dominated outer regi ons left over by mergers 
(e.g. NGC 5266, see lMorganti et al J 199% . To account for the dif- 
ferent behavior of R^ xy in elliptical and spiral galaxies, we chose a 
piecewise power-law with different powers for the two populations, 



iog(<H y ) = 



cf + uf T 



+ u 



if T < r* 
if T > 77 



' phy.l 



(ID 



where cf, uf, cf, uf are considered as the free parameters to be 
fitted to the data, and T* is at the intersection of the two regressions, 



..el 



+ uf 77 



' + uf Tp thus ensuring that R^^? remains a 
continuous function of T at T = T*. 

Another correlation exists between R g,l] '" iy and the total cold 

mol 

gas mass M gas or between R^ xy and the blue magnitude Mb- In 
fact, these two correlations are closely related due to the mutual 
correlation between M m and M B , and hence we shall restrict our 
considerations to the correlation between R^ xy and M g!ls . Accord- 
ing to the roughly monotonic trend visible in Fig. [5] we choose a 
power-law between R^ xy 



q 2 + k 2 log 



10 9 /r 2 M G 



phy.i 



(12) 



where q 2 , k 2 are free parameters. A somew hat similar depend ence 
was recently found between 7?^ xy and M m (Ker es et alj2003l) , but 



Understanding the H2/HI Ratio in Galaxies 7 



Model log(fi; 



;;il;ixy ^ 



-0.58+ 



-0.51 + 



k 
T* 



+ 0.18!g;g 
+0.12^ 

-0 14+ 010 
V - ^-0.07 

-0 12 +0 - 01 



-1.3+ 



-0.24+ 



-0 01+ - 25 

uul -0 .16 

+0 13+ 007 

-0 02+ 010 
u - -0.09 

-0 13 +002 

u - -0.02 

n 1 s+0.06 
-U.18_ 007 

-0 1 +L2 



O" data.i 
°" phy.i 



0.71 

0.39 



0.66 
0.27 



0.67 
0.30 



0.62 
0.15 



Table 4. The upper panel lists the most likely parameters and 1-cr confi- 



dence intervals of the four models R' 



galaxy 



a = o, 



. ,, 3). The bottom panel 

shows the rms-deviations <x data- j of the data from the model predictions and 
the estimated physical scatter o" p h y- i for each model i. 



this result is less conclusive, since R^^ y and M m are naturally cor- 
related by the definition of R^^, even if M ra and M H2 are com- 
pletely uncorrected. 

Finally, we shall introduce a fourth model (model 3) for ^' ',' xy 
that simultaneously depends on galaxy Hubble type and cold gas 
mass, 



log(iO 



+h log 



[ T (ifT<T;) 
(ifT>T;) 



(13) 



10 9 h- 2 M, 



+ °~ phy.3 , 



where cf, 
of T at T = 



c*, uf, kj, are free parameters and is defined as 
■? + w! p 7"j, thus making R^Y a continuous function 
'3. Comparing this model with models 1 and 2, will 
also allow us to study a possible degeneracy between model 1 and 
model 2 caused by a dependence between cold gas mass and galaxy 
Hubble type. 

The free parameters of the above models were determined by 
minimizing the rms-deviation between the model predictions and 
the 245 observed values of log(2? gal J xy ) (Appendix [A). Optimiza- 
tion in log-space is the most sensible choice since R^^ is subject 
to Gaussian scatter in log-space as will be shown in the Section 
14.31 The most probable values of all parameters are shown in Table 
[4] together with the corresponding 1-cr confidence intervals. The 
latter were obtained using a bootstrapping method that uses 10 4 
random half-sized subsamples of the full data set and determines 
the model-parameters for every one of them. The resulting distri- 
bution of values for each free parameter was approximated by a 
Gaussian distribution and its standard deviation cr was divided by 
V2 in order to find the 1-cr confidence intervals for the full data set. 
Note that in some cases the parameter uncertainties are coupled, 
i.e. a change in one parameter can be accommodated by changing 
the others, such that the model remains nearly identical. For models 
1 and 2, the best fits are displayed in Figs.[4]and[5]as solid lines. 

Table [4] also shows different scatters that will be explained in 
Sectionl4~3l 



4.3 Scatter and uncertainty 



The empirical values of ^^" y scatter around the model predictions 
according to the distributions shown in Fig. [6] (dashed lines). The 
close similarity of these distributions to Gaussian distributions in 
log-space (solid lines) allows us to consider the rms-deviations of 
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Figure 6. Distributions of the deviations between the observed values of 



log(S^ xy ) and the model-values log(fi^" J ) = 1, 3). Data points and 
dashed lines represent the actual distribution of the data; solid lines repre- 
sent Gaussian distributions with equal standard deviations. 



.galaxy. 



the data c data as the standard deviations of Gaussian distributions. 
This exhibits the advantage that cr dala can be decomposed in model- 
independent observational scatter cr obs and model-dependent phys- 
ical scatter cr n; 
( = 0,...,3. 



1 p hy via the square-sum relation cr 2 ata . = < r 



obs ^ ^~phy.i 



The major contribution to cr da , ai comes from observational 
scatter, as suggested by the close similarity of the different values 
of cr datai . Indeed, the observational scatter inferred from the R^T y 
values of the 22 repeated sources in our data is cr obs as 0.6. This 
scatter is a combination of CO-flux measurement uncertainties, un- 
certain CO/H2-conversions and Hi-flux uncertainties (in decreas- 
ing significance). Since c obs is only marginally smaller than cr dataii 
for all models, the estimation of the physical scatters cr phy ; (given 
in Table f4j is uncertain. Nevertheless, we shall include these best 
guesses of the physical scatter, when constructing the H2-MFS in 
Section|4T4l 



4.4 Recovering the H 2 -MF and model evidence 

Given a model for R^ xy ', H 2 -masses of arbitrary Hi-galaxies can 
be estimated. We shall apply this technique to the 4315 sources in 



the HIPASS catalog using our four models of R 



lalaxv 



0.....3. 



For each model, the resulting H 2 -catalog with 4315 objects will be 
converted into a H 2 -MF, which can be compared to our reference 
H 2 -MF derived directly from the CO-LF (Section|5). 



For the models R g ^{(T) and R^(M gss , T) Hubble types T 
were drawn from the HyperLeda database for each galaxy in the 
HIPASS catalog by means of the gal axy identifiers giv en in the op- 
tical cross- match catalog HOPCAT jDovle et alj|2005h . H 2 -masses 
were then computed via Mn 2 = R^JIT Mm, i = 0, 3. This equa- 
tion is implicit in case of the mass-dependent models f^ y (M gas ) 

and R^(.M m , T), where M gas = (M m + M H2 )/j3. All four models 
were applied with scatter, randomly drawn from a Gaussian distri- 
bution with the model-specific scatter cr phy ;, listed in Table|4] 

In order to reco nstruct a H 2 -MF for each model, we employed 
the 1/Vmax method dSchmidt|[l968l) , where V max was calculated 
from the analytic completeness function for HIPASS that depends 
on the HI peak flux density S u , the integrate d HI line flux S mt , and 
the flux limit of the survey jZwaan et alj2004h . After ensuring that 
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Figure 7. H2-MFS constructed from the HIPASS Hi-catalog using the dif- 
ferent phenomenological models for the HI/H2 ratio. The black dashed line 
is the reference H2-MF derived from the CO-LF in Section[5] 

we ca n accurately reproduce the HI-MF derived by IZwaan et al.l 
(2005), we evaluated the four H 2 -MFs (one for each model ^ m ^|* y ) 
displayed in Fig.[7](dots). The uncertainties of log(0H 2 ) vary around 
cr = 0.03 - 0.1. Each function was fitted by a Schechter function 
by minimizing the weighted rms-deviation (colored solid lines). 

The comparison of these four H2-MFS with the reference H2- 
MF derived from the CO-LF allows us to qualify the different mod- 



els R 



.galaxy 



0, 3, against each other. We ask: "What are the 

odds of model R^J^? against model R^* y if the reference H2-MF 
derived from the CO-LF is correct?" This question takes us back 
to the Bayesian framework of model selection applied in Section 
12.21 If the models are a priori equally likely, their odds are equal to 
the Bayes factor, defined as the ratio between the model evidences. 
When computing these evidences, we take the "observational" data 
d to be the reference H2-MF (with scatter), while the "model" data 
is the H 2 -MF reproduced by applying a model ^^" y to the HIPASS 
(vector) are listed in Table [4] for each 
r or model #ET, y ). The P r i° r probabil- 
ity density n(8\Mi) in the evidence integral of Eq. © is taken as 
the multi-dimensional parameter probability distribution function 
obtained from the 245 galaxies studied in Section POl (see Table 
@). The second piece in the evidence integral, i.e. the probability 
density p(d\9, Mi), is calculated as the product, 



data. The free parameters 
model (e.g. 



p(.d\0,Mi) = V\ = exp 



« - <C dcU ) 2 



2o- 2 



(14) 



where k labels the different bins of the H2-MF (Fig. [7), and tp^ 
and 0™ odc1 ' 1 respectively denote the differential mass densities of 
the reference H 2 -MF and the fL-MFs reconstructed from HIPASS 
i = 0, ...,3. cr denotes the combined 

imodel.i 



using the models R^y^, 



statistical uncertainties of <p^ f and <p™° aa "\ theoretically given by 
cr 2 = cr™ (2 + -™ odc1 - 1 _ However, we shall neglect the contribution of 
cr™ * 1 ' 1 , since cr£ ef is about 3-4 times larger due to the small size of 
the FCRAO sample of CO-galaxies compared to the HIPASS sam- 
ple of Hl-galaxies. Furthermore, we assume that a is independent 
of the bin k and adopt an average uncertainty equal to cr = 0. 15 dex. 
This is the mean scatter of the binned data of the reference H 2 -MF 
(see Fig. [3}. Assuming a constant scatter for the whole reference 
MF artificially increases the weight of the low and high mass ends, 
where the scatter is indeed closer to 0.3 dex, and reduces the weight 



Model 


pgalaxy 
mol.O 


,,«alaxv 
/< . . 
mol, 1 


D »alaxv 
K . , 
mo 1,2 


^galaxy 
mol, 3 


Nb. of free param. 
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lnfi 


0.0 


7.3 


8.2 


22 



Table 5. Comparison of different models for H2/HI-mass ratios of entire 
galaxies; the first row shows the number of free parameters, while the sec- 
ond row shows the model evidence in terms of the Bayes factor between 
that model and the constant model R gl 



ogalaxy 
mol,0 ' 



of the central part, where the scatter equals 0. 1 dex. We argue that 
this is a reasonable choice, since the central part of the reference 
H2-MF suffers most from systematical uncertainties of the X-factor 
and the low and high mass ends encode much of the physics that 
could discriminate our models for R^ xy against each other. In any 
case, the outcome of this evidence analysis is only weakly affected 
by the choice of scatter. 

The integration of the evidence integral is computationally ex- 
pensive: for each choice of model-parameters the following three 
steps need to be performed: (i) evaluation of the H 2 -masses for 
each galaxies in the HIPASS sample, (ii) computation of the H2- 
MF from that sample, (iii) computation of the product in Eq. t!4t . 
We applied a Monte Carlo method to sample the parameter spaces 
of the different models. About 10 6 integration steps had to be per- 
formed in total to reach a 2 per cent convergence of the Bayes fac- 
tors. 



The Bayes factor between each model R\ 



galaxy 



0,...,3, and 



^moLo' ' s shown in Table[5] We find strong evidence for all variable 
models (R^f, R^Jrl, ^maT) a g a inst the constant one (R^^), 
and there is even stronger evidence of the bilinear model (f?5^ y ) 
against all others. The H2-MF associated with this model is indeed 
the only one providing a simultaneous fit to the low and high mass 
ends of the reference MF (see Fig. |7J, and the agreement is very 
good (reduced^ 2 = 0.8). 

On a physical level, there are good reasons for the partial fail- 
ure of the other models in reproducing the extremities of the ref- 
erence H2-MF. Model R^7(T) overestimates the space density of 
galaxies with high H 2 -masses by overestimating R^^ for the gas- 
richest early-type spiral galaxies. In reality, the latter have a very 
low molecular fraction (see data, model R^J^ , theory in Section|5]l, 
but they are a minority within otherwise gas-poor but molecule-rich 
early-type spirals. Hence, a model depending on Hubble type alone 
is likely to miss out such objects, resulting in an increased density 
of high H 2 -masses. While model R^^(M gas ) overcomes this is- 
sue and produces the right density of high H2-masses, it fails by 
a factor 3 - 4 in the low-mass end (Mu 2 < 10 8 M o ). This is a di- 
rect manifestation of assigning high molecular fractions to all gas- 
poor galaxies, which neglects small young spirals with a dominant 
atomic phase. Finally, model R^^ seems to suffer from limitations 
at both ends of the H 2 -MF 

The clear statistical evidence for model 3 shall be supported 
by the theoretical derivation of /? gdl ' lxy presented in Section[5] 



5 THEORETICAL MODEL FOR THE H 2 /HI-MASS 
RATIO 

So far, we have approached the galactic H 2 /HI-mass ratios 
^moi* 5 ' with a set of phenomenological models, limited to the 
local Universe. By contrast, we have recently derived a phys- 
ical model for the H2/HI-ratios i n regular galaxies, whic h po- 
tentially extends to high redshift dObreschkow et all [2009), This 
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model relies on the theoretically and empirically established re- 
latio n between interstellar gas pressure and local molecular frac- 
tion dElmegreerJI 19931: )Blitz & Rosolowskvll2006l : lKrumholz et al.l 
l2009l : lLerov et al.l2008l) . In this section, we will show that the phys- 
ical model predicts H 2 /HI-ratios consistent with our phenomeno- 
logical model 3 given in Eq. l |13t . Hence, the physical (or "theo- 
retical") model provides a reliable explanation for the global phe- 
nomenology of the H 2 /HI-ratio in galaxies. 

5.1 Background: the i? rao |-pressure relation 

Understanding the obs erved continuous v ariation of R mo \ within in- 
dividual galaxies (e.g. iLerov et al.ll2008l) requires some explana- 
tion, since, fundamentally, there is no mixed thermodynamic equi- 
librium of HI and H2. To first order, the ISM outside molecular 
clouds is atomic, while a cloud-region in local thermodynamic 
equilibrium (LTE) is either fully atomic or fully molecular, depend- 
ing on the local state variables. The apparent continuous variation 
of R mol is the combined result of (i) a non-resolved conglomera- 
tion of fully atomic and fully molecular clouds, (ii) clouds with 
molecular cores and atomic shells in different LTE, and (iii) some 
cloud regions off LTE with actual transient mixtures of HI and 
H2. However, a time-d ependent model for off-equilibrium clouds 
(Goldsmith et al. 2007) revealed that the characteristic time taken 
between the onset of cloud compression and full molecularization 
is of the order of 10 7 yrs, much smaller than the typical age of 
molecular clouds, and hence the fraction of these clouds is small. 
Therefore, averaged over galactic parts (hundreds or thousands of 
clouds), R mo \ is dictated by clouds in LTE, entirely defined by a 
number of state variables. 

A theoretical fra me explo i ting th e LTE of molecular clouds 
was introduced by lElmegreerJ l|l993[), who considered an ideal- 
ized double population of homogeneous diffuse clouds and isother- 
mal self-gravitating clouds, both of which can have atomic and 
molecular shells. In this model the molecular mass fraction f mol = 
dMH 2 /d(Mni + Mn 2 ) of each cloud depends on the density profile 
and the photodissociative radiation density from stars j, corrected 
for self-shielding by the considered cloud, mutual shielding among 
different clouds, and dust extinction. Since the sh i eldin g from this 
radiation depends on the gas pressure, Elmegreen d 19931) finds that 
/moi essentially scales with the external pressure P and photodisso- 
ciative radiation density j, approximately following f mo \ oc P 22 j' 1 
with an asymptotic flattening towards f mol = 1 at high P and low 
j. This implies approximately R mo \ = dMn 2 /dMm K P j ■ As- 
suming that j is proportional to the surface density of stars E staIS 
and that the ste llar velocity dispersion cr stars vari e s radi a lly as 
IWong & Blitz! J2002h and iBlitz & Rosolowskvl d2004l 120061) find 
roughly j oc P and h ence R m0 \ oc P a with a = 1.2. Recently, 
iKrumholz et al.l (2009) have presented a more elaborate theory con- 
cluding that a a 0.8. However, the exponent a remains uncertain, 
thus requiring an emp irical determination. 

Observationally, IBlitz & Rosolowskvl d2004 l200rj) were the 
first ones to reveal a surprisingly tight power-law relation between 
pressure and molecular fraction based on a sample of 14 nearby 
galaxies including dwarf galaxies, Hi-rich galaxies, and H 2 -rich 
galaxies. Per haps the richest ob servational study published so far 
is the one bv lLerov et al.l d2008I) . wh o analyzed 23 gala xies of The 
HI Nearby Galaxy Survey (THINGS. IWalter et alj200"j) , for which 
H 2 -densities had been derived from CO-data and star formation 
densities. This analysis confirmed the power-law relation 

Rmol = {P/P.y, (15) 



where P is the local, kinematic midplane pressure of the gas, and 
P, and a are free parameters, whose best fit to the data is given by 
P. = 2.35 • 10 13 Pa and a = 0.8. 



5.2 Physical model for the H 2 /HI-ratio in galaxies 

We shall now consider the consequence of the model given in 
Eq. J 1 5b for the H2/HI-ratio of entire galaxies. To this end, we adopt 
the models and methods presented in lObreschkow et al. I J2009\ re- 
stricting this paragraph to an overview. 

First, we note that most cold gas of regular galaxies is nor- 
mally contained in a disc. This even applies to bulge-dominated 
early-type galaxies, such as suggested by recently pr esented CO- 
maps of five nearby elliptical galaxies (Young 2002). Hence, the 
HI- and H 2 -distributions of all regular galaxies can be well de- 
scribed by surface density profiles Ehi('') and EjfeW- We assume 
that the disc is composed of axially symmetric, thin layers of stars 
and gas, which follow an exponential density profile with a generic 
scale length Jdj SC , i.e. 

E^ c s (r) ~ E gas (r) ~ E H ,(r) + Sh 2 W ~ exp(-r/r dlsc ), (16) 

where r is the galactocentric radius and £ denotes the mass col- 
umn densities of the different components. Next, we adopt the phe- 
nomenological relation of Eq. l |15t , i.e. 

= [P(r)/P,] a , (17) 

and s ubstitute the kinematic midplane pressure P(r) for ( Elmegreen 

Il989h 

P(r) = n -G X gas (r)(s gas (r) + /l£S(r)), (18) 

where G is the gravitational constant and / = a" gas ,z/0" starSiZ is the 
ratio between the vertical velocity dispersion s of gas and stars. We 
adopt / = 0.4 according to lElmegreenl ( 1989). 

Eqs. l|16| [TTJ can be solved for Zhi(>\) an d £h 2 (')- In 
lObreschkow et al.ld2009h . we demonstrate that the resulting surface 
profiles are consistent with the empiric al data of the two nearby spi- 
ral galaxies NGC 5055 and NGC 5 194 dLerov et alj2008r) . Integrat- 
ing ShiM an d 2h, ( r ) over the exponential disc gives the gas masses 
Mm and Mh 2 , hence providing an estimate of their ratio #jL ] xy . An- 
alytically, R^* y is given by an intricate expression, which is well 
approximated (relative error < 0.05 for all galaxies) by the double 
power-law 

= ( 3 - 44 ^°- 5 ° 6 + 4.82iCr 1 -° 5 T 1 - ( 19 ) 

where 

K™, = [ 1 1 .3 m 4 kg- 2 r di t. M gas (M gas + 0.4 M^ s )f 8 . (20) 

is a dimensionless parameter, which can be interpreted as the 
H 2 /HI-ratio at the center of a pure disc galaxy. For typical cold gas 
masses of average galaxies (M gas = 10 8 - 10'°M o ) and correspond- 
ing stellar masses and scale radii, R?j calculated from Eq. J20t 
varies roughly between 0. 1 and 50. Hence, ^^„" y given in Eq. ([19} 
varies roughly between 0.01 and 1. 

In summary, Eqs. U9I I20t represent a theoretical model for 
R^ xy , which uses three input parameters: the disc stellar mass 
M f!^, the cold gas mass M ga s , and the exponential scale radius r disc 
(see lobreschkow et alj2009l for a detailed discussion). 
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Figure 8. Visualization of the function ^^"oi^^mol^ T ne solid line repre- 
sents the nearly exact function given in Eq. 1191 . while the dashed line is 
the power-law fit of Eq. i2i\ . 



5.3 Mapping between theoretical and phenomenological 
models 

We shall now show that our theoretical model for galactic H2/HI- 
mass ratios given in Eqs. d 1 91 l20t closely matches the best phe- 
nomenological model given in Eq. dl3t . The mapping between the 
two models uses a list of empirical relations derived from obser- 
vations of nearby spiral galaxies, and hence the comparison of the 
models is a priori restricted to spiral galaxies in the local Universe. 

First, we note that Eq. d 1 91 > can be well approximated by the 
power-law 

O-O.l^o, 08 . (2D 

As shown in Fig. [8] this approximation is accurate to about 10 per 
cent over the whole range R?. = 0.1, ...,50, covering most regular 
galaxies in the local Universe. 

Substituting R9. in Eq. i2U for Eq. l |20| l, yields the approxi- 
mate relation 

nO.64 



/? gal " y = [0.31m 4 kg- 

mol,th [ B 



, (Mga, 



iff 



(22) 



In order to compare the theoretical model of ^ gjl ' lxy to the em- 



pirical models of Section 1421 we need to eliminate the formal de- 
pendence of RT^*1 on r «iisc an d ^tars- To this end, we use two ap- 
proximate empirical relations, derived from samples of nearby spi- 
ral galaxies (see Appendix lBl>. 



log I 



log I 



J " stars 

/r 2 M Q 

Misc 

/r'kpc 



7i + on log 



72 + a% log 



2 • 10" hr 2 M, 



5 • 10" h- 2 M, 



+ 6T, 



(23) 



(24) 



where T = (10— 7^/16 is the normalized Hubble type, which varies 
between T = (pure disc galaxies) to T = 1 (pure spheroids). 

The parameters corresponding to the best x 1 nt (Appendix 
HJl are a, = 1.46 ± 0.1, y, = 9.80 ± 0.05, a 2 = 0.45 ± 0.05, 
y 2 = 0.97 ± 0.05, 6 = -1.07 ± 0.1. The given intervals are the 1- 
cr confidence intervals of the parameters; they do not characterize 
the scatter of the data. The units on the right hand side of Eqs. d23l 
124b were chosen such as to minimize the correlations between the 
uncertainties of a, and y\. 

Physical reasons for the empirical relations in Eqs. J23II24I ) are 
discussed in Appendix [B] Substituting Eqs. J23II24I ) into Eq. d22t 
reduces irj? to a pure function of M gas and T of the form 

log[CS(^s,r)] = log[<-(M g a S ,10)] (25) 




log(M gas /[10 9 /,- 2 M o ]) 

Figure 9. Relation between R®£% and M gas for flat disks (T = 10). The 
solid line represents the relation obtained from Eq. (22), when expressing 
''disc and AfiS? as functions of M gas using Eqs. 12311241 . The shaded zone 
represents the 1 -cr uncertainty implied by the uncertainties of the empirical 
parameters in Eqs. j23H24t . The dashed line represents the best power-law 
fit for the displayed mass interval as given in Eq. i26\ . 



+5(0.167-1.6), 



where R^hiM^, 10) is the theoretical Hj/HI-ratio of a pure disc 
galaxy, i.e. T = 10. The function R^^l(M gas , 10) is displayed in 
Fig.|9]together with the 1-cr uncertainty implied by the uncertainties 
of the four parameters ct\, a%, 71, 72- We approximate this relation 
by the power-law 



10) =c+ J -log 



10 9 hr 2 M e 



(26) 



The parameters minimizing the rms-deviation on the mass-interval 
log(M gas /[M /T 2 ]) = 7.5 - 10.5 are c = -1.79 ± 0.04 and s = 
-0.24 ± 0.05. The given uncertainties approximate the propagated 
uncertainties of a u a 2 , 71, 72. 

The simplified theoretical model for the lU/HI-ratio given in 
Eqs. l !25ll26t exhibits exactly the formal structure of our best phe- 
nomenological model 3. Setting R ^ i^(Mgas> T) in Eq. l |25t equal to 
/?^j y (M gas , T) in Eq. J 1 3t for spiral galaxies, yields the following 
mapping between the theoretical and empirical model-parameters, 



c- 1.65, 
0.165. 



(27) 



The probability distributions of the empirical model- 
parameters on the left hand side of Eqs. J27t were derived in 
Section [4] and their l-cr uncertainties are given in Table [4] The 
corresponding probability distributions of the theoretical model- 
parameters on the right hand side of Eqs. i27l can be estimated 
from the Gaussian uncertainties given for the parameters c, s, 6. 
The empirical and theoretical parameter distributions are compared 
in Fig.[l0]and reveal a surprising consistency. 



6 DISCUSSION 

6.1 Theoretical versus phenomenological model 

The dependence of ^j^[ Xy on galaxy Hubble type T and cold gas 
mass Mgas was first considered on a purely phenomenological level, 
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Figure 10. Probability distributions of the three parameters in our model 
3 (Eq. 113) for the H2/HI-mass ratio fij^ xy of spiral galaxies. Solid lines 
represent phenomenologically determined probability distributions given in 
Table [4] dashed lines represent the corresponding theoretical probability 
distributions, obtained when using Eqs. (27) with the respective distribu- 
tions for c, and 6. 



and described by the empirical models in Section [4] The best em- 
pirical model for spiral galaxies could be quantitatively reproduced 
by the subsequently derived theoretical model for regular galaxies 
in Section[5] Hence, the latter provides a tool for understanding the 
variations of R^ xy ■ 

In fact, according to Eq. {22}, R^ xy seems most directly dic- 
tated by the scale radius rdi sc and the masses M gas and M^ s . The 
dependence of R^ xy on T is clearly due to the trend for smaller 
values of rdi sc (for a given mass) in bulge-rich galaxies. Several 
physical reasons for the influence of the bulge on are men- 
tioned in Appendix |B2] 

From Eq. d22b . one might naively expect that ^f^ xy and M gag 
are positively correlated. However, the disc scale radius rdi sc in- 
creases with M g!ls as oc M^ 6 by virtue of Eqs. J23II24I ). Tak- 
ing this scaling into account, the H2/HI-ratio R^ xy effectively de- 
creases with increasing M m . The physical picture is that more mas- 
sive galaxies are less dense due to their larger sizes, and hence their 
molecular fraction is lower. 

The 'best' phenomenological model is by definition the one 
that, when applied to the galaxies in the HIPASS sample, exhibits 
the H 2 -MF that best fits the reference H 2 -MF derived from the CO- 
LE The close agreement between the best model defined in this 
way and the theoretical model therefore supports the accuracy of 
the CO-LF iKere s etalj|2003h . which could a priori be affected 
by the poo rly characterized c ompleteness of the CO-sample. Con- 
firminglv. lKeres et al.l (2003) argued that the CO-LF does not sub- 
stantially suffer from incompleteness by analyzing the FIR-LF pro- 
duced from the same sample. 



6.2 Brief word on cosmic evolution 

galaxy 



The theoretical model given in Eqs. d 1 91 |20t potentially ex- 

tends to high redshift, as it only premises the invariance of the re- 
lation between pressure and R mo ] and a few assu mptions with weak 
depen dence on redshift (but see discussion in Obre schkow et al.l 
2009). However, we emphasize that the transition from the theoret- 



ical model /fjj^jj to the phenomenological model R^. , 
of relations extracted from observations in the local Universe. Most 
probably R^7 underestimates the molecular fraction at higher 
redshift, predominantly due to the evolution in the mass-diameter 



^galaxy 



uses a set 



relation of Eq. d24t . Indeed, scale radii are smaller at higher red- 
shift for iden tical masses, thus incre asing the pressure and molec- 
ular fraction. Bou wens et all d2004l) found r disc oc (1 + z) 1 from 
observations in th e Ultra Deep F ield, consistent with the theoret- 
ical prediction by IMo et alj Jl998t) . According to Eq. d22t . where 



^galaxy ^ 2 . 6 



this implies R^ xy oc (1 + z) 2b - In other words, 



mol " ' disc 

the phenomenological model 3 (Eq. 1131 ) for the H2/HI-mass ratio 
should be multiplied by roughly a factor (1 + z) 2 6 . However, this 
conclusion only applies if we consider galaxies with constant stellar 
and gas masses. For the cosmic evolution of the universal H2/HI- 
ratio ^™ vclsc , we also require a model for the evolution of the stellar 
and gas mass functions, and it may even be important to consider 
different scenarios for the evolution of the scale radius for different 
masses. A more elaborate model for the evolution of R l Z" crx 



be ob tained from cosmological simulations (e.g. lObreschkow et al.l 
2009 and forthcoming publications). 



6.3 Application: The local cold gas-MF 

We finally apply our best phenomenological model for the H2/HI- 
mass ratio (i.e. R^*J given in Eq. 1131 ) to derive an integral cold 
gas-MF (HI+H 2 +He) from the HIPASS catalog. In fact, the cold 
gas-M F cannot be inferred solely from the HI-MF (e.g. Zwa an et"al] 
2005) and the H 2 -MF (e.g. Section [3}, but only from a sample of 
galaxies with simultaneous HI- and H2-data. Presently, there is no 
such sample with a large number of galaxies and an accurate com- 
pleteness function. Therefore, we prefer using the HIPASS data, 
which have both suffic ient size (4 315 galaxies) and well described 
completeness dZwaan et alj2004h . and we estimate the correspond- 
ing H 2 -masses using our model ■ Details of the computation 
of the H2-masses were given in Section l4~4l 

The resulti ng cold gas-MF is s hown in Fig. QT] together with 
the HI-MF from lZwaan et alj J2005I) and the reference H 2 -MF de- 
rived in Section[3] The displayed continuous functions are best fit- 
ting Schechter functions. The respective Schechter function param- 



eters for the cold gas-MF are M* = 7.21 • 10 9 /T 2 M e 



•1.37, 



and (f>* = 0.01 14 h 3 Mpc 3 . The total cold gas density in the 
local Universe derived by integrating this Schechter function is 
Q gas = 4.2- 10~ 4 h~ x , closely matching the value (4.4+ 0.8)- 10'" hr 1 
obtain ed when summing up the empirical Hi-density (Zwa an et alj 
2005), the H2-density (Section [3}, and the corresponding He- 
density. 



7 CONCLUSION 

In this paper, we established a coherent picture of the H 2 /HI-ratio 
in galaxies based on a variety of extragalactic observations and the- 
oretical considerations. Some important jigsaw pieces are: 

(i) Measurements of the Jf -factor (summarized in Bose lli et al.1 
l2002h were combined with more recent CO-flux measurements and 
extinction-corrected optical M B -magnitudes, resulting in a working 
model for X. 

(ii) This model for X was applied to the CO-LF bv lKeres et al.1 
(2003) in order to derive the first local H 2 -MF based on a variable 
X-factor. 

(iii) Nine samples of local galaxies (245 objects in total) with 
simultaneous measurements of Mm and Leo were combined to fit 
a set of empirical models for galactic H 2 /HI-mass ratios R^ xy ■ 

(iv) These models were applied to the large Hl-sample of the 
HIPASS catalog, which permitted the derivation of a H2-MF for 
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Figure 11. Filled squares represent the integral cold gas-MF (HI+tb+He) 
derived from the HIPASS data using our best phenomenological model for 
the H2/HI-mass ratio ( Eq. 1131 ; empty squares represent the observed HI- 
MF (Zwaan et al. 2005) and empty circles represent our best estimate of the 
H 2 -MF (Section[3J. Solid lines are best fitting Schechter functions. 

each model for R^'^ ■ A comparison of these H 2 -MFs with the one 
derived directly from the CO-LF allowed us to determine the sta- 
tistical evidence of each model and to uncover a clear 'best model'. 

(v) Based on the relation between pressure and the local H2/HI- 

1 ■ 11 — ■ — 1 

ratio ^ mol dLerov et al . 2008), we established a theoretical model 

for the H2/HI-ratio R^™ y of regular galaxies, which potentially ex- 
tends to high redshifts. 

(vi) We could show that the best empirical model for ^j^ xy 
found before is an excellent approximation of the theoretical model 
in the local Universe. 

The factual results standing out of this analysis are 

(i) an empirica l H 2 -MF obtained by combining the CO-LF of 
iKeres et al. (2003]) with a variable Z-factor (see Fig. [3] and param- 
eters in Table [3}, 

(ii) an empirical model for R^*^ (Eq.ll3t. which accurately re- 
produces the above H2-MF, when applied to the Hl-sample of the 
HIPASS catalog, 

(iii) a theoretical model for R^T y (Eqs. |19l|20"l >. which provides 
a source for physical understanding and generalizes to high red- 
shift, 

(iv) a quasi-empirical integral cold gas-MF (HI+H2+He) based 
on the HIPASS data. 

Self-consistency argues in favour of the interlinked picture es- 
tablished in this paper. However, all quantitative results remain sub- 
jected to the uncertainties of the X-factor. The latter appears as a 
scaling factor, affecting in the same way the reference H 2 -MF de- 
rived from the CO-LF, the phenomenological models of R^ xy and 
hence the H2-MFS derived from HIPASS, as well as the P-R m0 \ re- 
lation and thus the theoretical model for R^ xy ■ In the future it may 
therefore be necessary to re-scale the quantitative results of this pa- 
per using a more accurate determination of X. 
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1989). These 267 objects contained 22 repeated galaxies. In each 
case of repetition, the older reference was removed, such as to re- 
main with the 245 distinct sources listed in Table IA1I The CO- 
luminosities were homogenized by making them independent of 
different .X-factors and Hubble constants. All other properties listed 
in the table were taken from homogenized reference catalogs, such 
as described in Section RTl 



APPENDIX A: HOMOGENIZED DATA 

This section presents the data (245 galaxies) used for the derivation 
of the models of R^* y in section [4] 

CO-luminosities were drawn from 10 small er samples: 17 



nearby (< 10 Mpc ) lenti culars and ellipticals (Welch & Sage 



20031: ISage & W elch 2006), 4 late-type spirals dMatthews et al 



20051). 68 isolated late-ty pe spirals dSautv et ai1l2003l). 6 ellipticals 



(Georgaka kis et al]|200ll) . 17 spirals of all types dAndreani et alj 
1 19951) , 48 near by (< 10 Mpc) spirals of all types dSagelll993h 



12 ellipticals i Lees et al.1 Il991), 18 lenticula rs and ellipticals 
dThronson etauT l989), 77 spirals of all types dYoung & Kn ezek 



D. Obreschkow and S. Rawlings 



Table Al: Homogenized galaxy sample based on data drawn from the literature. 
T is the numerical Hubble type (see online help of the HyperLeda database), 
D[ the luminosity distance, M B is the extinction corrected absolute blue mag- 
nitude, and X is the variable X-factor derived from M B (Eq. [8) without addi- 



J2003b: Saee & Welch 


( 


2006 


), mlMatthews et alJd2005b. mlSautv 


etal. 12003), 


[4] Georgakakis et al. 




2001 


). T51 Andreani et al. 1 1995b. T61 Sage 


: H993), [7] 


ILees et all ill 99 lb. T81 


Thronson et alJdl989h. T9llYoune & Knezek 


(1989). 



Object 


T 


D, //r'Mpc 


M B - 5 log h 


X 


log(M H2 /X/7- 2 M Q ) 


Ref. H 2 


log(M ffl //T 


NGC 404 


-2.8 


1.7 


-15.86 


6.66 


6.06 


1 


7.51 


NGC 2787 


-1.1 


9.5 


-18.87 


2.21 


6.58 


! 


8.58 


NGC 3115 


-2.8 


6.4 


-19.27 


1.91 


5.60 


! 


6.75 


NGC 3384 


-2.7 


9.2 


-19.06 


2.06 


5.87 


1 


5.94 


NGC 3489 


-1.3 


7.7 


-18.45 


2.58 


6.12 


j 


6.46 


NGC 3607 


-3.1 


10.3 


-19.23 


1.94 


8.34 


! 


6.93 


NGC 3870 


-2.0 


9.9 


-16.56 


5.15 


7.44 




8.08 


NGC 3941 


-2.0 


11.0 


-19.04 


2.08 


7.15 


! 


8.81 


NGC 4026 


-1.8 


12.1 


-18.82 


2.25 


7.27 


1 


7.86 


NGC 4150 


-2.1 


6.8 


-17.66 


3.44 


6.91 


j 


6.88 


NGC 4203 


-2.7 


12.7 


-18.86 


2.22 


6.21 


! 


8.41 


NGC 4310 


-1.0 


10.8 


-16.86 


4.61 


6.96 


! 


7.10 


NGC 4460 


-0.9 


7.3 


-17.04 


4.32 


6.45 




8.26 


NGC 4880 


-1.5 


14.8 


-17.92 


3.13 


6.27 


1 


6.02 


NGC 7013 


0.5 


9.6 


-18.79 


2.28 


7.30 


j 


8.70 


NGC 7077 


-3.9 


12.0 


-16.13 


6.03 


6.09 


1 


7.60 


NGC 7457 


-2.6 


9.6 


-18.29 


2.73 


5.85 


! 


5.88 


NGC 100 


5.9 


9.0 


-17.61 


3.51 


5.91 


2 


8.87 


UGC 2082 


5.9 


7.7 


-17.72 


3.37 


5.89 


2 


8.80 


UGC 3137 


4.2 


12.5 


-17.05 


4.30 


6.20 


2 


9.11 


UGC 6667 


6.0 


12.1 


-17.06 


4.29 


5.73 


2 


8.54 


UGC 5 


3.9 


74.4 


-20.98 


1.02 


8.76 


3 


9.82 


NGC 7817 


4.1 


24.1 


-20.42 


1.25 


8.44 


3 


9.30 


IC 1551 


3.6 


136.0 


-22.17 


0.66 


8.94 


3 


9.34 


NGC 237 


4.5 


42.0 


-19.81 


1.57 


8.53 


3 


9.75 


NGC 575 


5.3 


32.3 


-19.10 


2.03 


8.04 


3 


9.18 


NGC 622 


3.4 


52.1 


-19.93 


1.50 


8.24 


3 


9.54 


UGC 1167 


5.9 


43.6 


-19.18 


1.97 


8.85 


3 


9.61 


UGC 1395 


3.1 


52.3 


-19.90 


1.51 


8.43 


3 


9.25 


UGC 1587 


3.7 


57.4 


-20.38 


1.27 


7.86 


3 


9.59 


UGC 1706 


5.8 


49.4 


-19.82 


1.56 


7.96 


3 


9.17 


IC 302 


4.1 


59.6 


-21.33 


0.90 


8.43 


3 


10.19 


IC 391 


4.9 


18.3 


-18.91 


2.18 


7.46 


3 


8.89 


UGC 3420 


3.1 


54.5 


-20.96 


1.03 


8.03 


3 


10.01 


UGC 3581 


5.2 


53.2 


-20.30 


1.31 


8.24 


3 


9.56 


NGC 2344 


4.4 


11.3 


-17.91 


3.14 


6.73 


3 


8.66 


UGC 3863 


1.1 


62.2 


-20.53 


1.20 


8.32 


3 


9.30 


UGC 4684 


7.2 


24.9 


-17.92 


3.13 


6.82 


3 


9.11 


NGC 2746 


1.1 


73.7 


-20.65 


1.15 


8.65 


3 


9.64 


UGC 4781 


5.9 


14.4 


-16.54 


5.19 


6.46 


3 


8.90 


UGC 5055 


3.1 


79.4 


-20.19 


1.36 


8.79 


3 


10.02 


NGC 2900 


5.9 


54.3 


-19.51 


1.75 


8.57 


3 


9.69 


NGC 2977 


3.2 


33.5 


-19.95 


1.49 


8.31 


3 


8.83 


NGC 3049 


2.5 


15.0 


-17.86 


3.20 


7.24 


3 


8.86 


IC 651 


8.2 


45.2 


-20.37 


1.28 


8.54 


3 


9.53 


NGC 3526 


5.2 


14.5 


-18.68 


2.37 


7.73 


3 


8.64 


UGC 6568 


8.2 


60.8 


-19.86 


1.54 


8.12 


3 


9.14 


UGC 6769 


3.0 


88.2 


-20.66 


1.15 


9.10 


3 


9.96 


UGC 6780 


6.4 


17.3 


-16.79 


4.73 


7.29 


3 


9.28 


UGC 6879 


7.1 


24.1 


-18.78 


2.28 


7.82 


3 


8.83 


UGC 6903 


5.9 


19.3 


-17.69 


3.40 


7.46 


3 


9.07 
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Object 


T 


A/[/r'M P c] 


M B - 5 log h 


X 


log(M H2 /[X/r 2 M Q ]) 


Ref. H 2 


log(M ffl /|7i : 


NGC 4348 


4.1 


20.3 


-19.49 


1.76 


8.10 


3 


9.01 


NGr 4617 


3.1 


49.6 


-20.70 


1.13 


8.56 


3 


9.90 


NGC 4635 


6.5 


10.9 


-17.28 


3.96 


6.73 


3 


8.23 


NGC 5377 


1.1 


20.6 


-19.83 


1.55 


7.81 


3 


8.91 




2.4 


76 

Z.U.U 


-19 54 


1 73 


7 60 


3 


9 74 


NGC 5584 


5.9 


17.1 


-19.06 


2.06 


7.22 


3 


9.27 


NGr 5690 


5.4 


18.4 


-19.88 


1.53 


8.15 


3 


9.33 


NGC 5768 

» > v J v .J / UO 


5.3 


20.3 


-18.74 


2.32 


7.90 


3 


9.11 


NGr 577? 


3.1 


52.3 


-20.41 


1.26 


8.25 


3 


9.49 


NGC 59 1 3 


1.3 


20.8 


-19.00 


2.1 1 


8.22 


3 


8.44 




1 9 
l.y 


20.1 


-19 00 


2.11 


7 73 


3 


9 76 


IC 1231 


5.8 


55.9 


-20.71 


1.13 


8.04 


3 


9.14 


TIGr 10699 

I V J V 1UU77 


4.4 


65.5 


-20.19 


1.36 


8.60 


3 


9.11 


TIGr 10743 


1.1 


27.2 


-18.75 


2.31 


7.52 


3 


8.78 


NGC 6347 


3.1 


64.3 


-20.46 


1.23 


8.57 


3 


9.48 


TIGr 10862 


5.3 


18.2 


-17.81 


3.26 


7.21 


3 


9.07 


NGr 6389 

. > V J V UJ07 


3.6 


33.1 


-20.37 


1.28 


8.30 


3 


9.93 


TIGr 1 1058 


3.2 


50.6 


-20.48 


1.22 


8.51 


3 


9.40 


NGC 6643 


5.2 


17.8 


-20.31 


1.30 


8.35 


3 


9.27 


NGr 67 1 1 

1 i V J V \J / 1 1 


4.0 


50.1 


-20.18 


1.37 


8.77 


3 


9.14 


TIGr 1 1 635 

I. . V J V 1 1UJJ 


3.7 


51.8 


-21.05 


0.99 


8.95 


3 


9.88 


TIGr 1 1773 


3.1 


50.1 


-19.87 


1.53 


8.48 


3 


9.57 


NGr 7056 


3.6 


55.8 


-20.53 


1.20 


8.67 


3 


9.11 


NGC 7 1 56 


5.9 


40.8 


-20.12 


1.40 


8.43 


3 


9.32 


T IGr 1 1 87 1 

l v J v 110/1 


3.1 


82.9 


-20.38 


1.27 


9.22 


3 


9.43 


NGr 7378 


2.1 


29.2 


-19.31 


1.88 


8.34 


3 


9.45 


NGr 7428 


1.1 


31.0 


-18.85 


2.23 


7.72 


3 


9.44 


TIGr 12304 

V V I V li.JU^ 


5.2 


35.3 


-19.40 


1.82 


8.01 


3 


8.88 


TIGr 12372 


4.0 


57.7 


-19.94 


1.49 


8.65 


3 


9.49 


NGC 7514 

1 *1 v I v / jit 


4.3 


51.1 


-20.62 


1.16 


8.18 


3 


9.16 


TIGr MAI A 


1.1 


53.5 


-20.53 


1.20 


8.80 


3 


8.87 


NGr 7664 

. N Vlv lUU^ 


5.1 


36.3 


-20.03 


1.44 


8.51 


3 


9.91 


TIGr 17646 


3.0 


83.7 


-20.84 


1.07 


8.68 


3 


9.70 


NGr 7712 


1.6 


31.9 


-18.94 


2.15 


7.84 


3 


9.10 


Tr 1 508 


7.2 


43.8 


-20.07 


1.42 


8.45 


3 


9.75 


TIGr 17776 

I . v J V. 1Z./ /U 


3.0 


51.8 


-19.88 


1.53 


8.31 


3 


9.99 


Tr 5355 

Iv, JJJJ 


5.7 


50.8 


-19.56 


1.72 


8.26 


3 


9.05 


TIGr 17840 


-1.8 


71.3 


-20.27 


1.32 


7.97 


3 


9.43 


NGr 7673 


2.0 


57.2 


-20.59 


1.18 


9.02 


4 


9.01 


NGr 7865 


-4.1 


26.0 


-20.01 


1.46 


7.35 


4 


8.79 


NGr 392 1 

- iviv j y 1 


0.0 


61.9 


-21.00 


1.01 


8.82 


4 


9.46 


NGr 4649 


-4.6 


12.1 


-20.70 


1.13 


7.15 


4 


8.35 


NGr 7252 


-2.1 


47.0 


-20.73 


1.12 


8.83 


4 


9.29 


NGr 7727 

1 1 V J V. / / — / 


1.1 


17.9 


-19.98 


1.47 


7.27 


4 


8.45 


NGr 142 

* > v IV 1 "-rz. 


3.1 


81.4 


-20.46 


1.23 


9.36 


5 


9.43 


IC 1553 


5.4 


28.0 


-18.75 


2.31 


7.69 


5 


9.10 


FSO 473-27 

ljkl\7 ^ 1 *J i-i 1 


4.4 


193.6 


-21.12 


0.97 


9.78 


5 


9.75 


NGC 232 

i > V I V Z.JZ 


1.1 


66.7 


-19.82 


1.56 


9.50 


5 


9.21 


FSO 475-16 


2.1 


70 7 


-20.44 


1.24 


9.01 


5 


9.74 


NGr 578 


5.0 


14.8 


-19.73 


1.61 


7.97 


5 


9.52 


J_.OvJ "4 / O-O 


zL 1 


S7 6 


70 6zL 


1 1 6 
1.10 


8 Q6 


< 

j 


Q 7^ 


NGC 1187 


5.0 


12.2 


-19.39 


1.83 


8.69 


5 


9.33 


NGr 1306 


2.8 


12.7 


-16.85 


4.63 


7.34 


5 


8.64 


NGC 1385 


5.9 


13.1 


-19.56 


1.72 


8.59 


5 


9.07 


ESO 549-23 


1.2 


40.8 


-19.42 


1.81 


8.45 


5 


8.88 


ESO 483-12 


0.3 


41.0 


-19.18 


1.97 


8.27 


5 


8.83 


NGC 1591 


1.9 


39.5 


-19.65 


1.66 


8.51 


5 


9.04 


NGC 7115 


3.4 


34.1 


-19.53 


1.73 


8.26 


5 


9.52 


NGC 7225 


-0.5 


47.9 


-20.09 


1.41 


9.29 


5 


9.07 


NGC 7314 


4.0 


13.2 


-19.71 


1.62 


8.05 


5 


9.24 
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Object 


T 


A/[/r'Mpc] 


M B - 5 log h 


X 


log(M H2 /[X/r 2 M Q ]) 


Ref. H 2 


log(M ffl /|7i : 


NGC 628 


5.2 


6.9 


-19.84 


1.55 


8.55 


6 


9.73 


NGC 67? 


6.0 


5.1 


-19.03 


2.08 


6.60 




9.07 


NGr 891 


3.0 


6.7 


-19.43 


1.80 


8.97 




9.72 


NGC 925 


7.0 


6.6 


-19.32 


1.87 


8.04 




9.57 


NGr 1058 


5.3 


6.3 


-17.78 


3.29 


7.42 




8.93 


NGC 1 560 


7.0 


2.3 


-15.91 


6.53 


5.88 




8.47 


NGr 7403 


6.0 


3.2 


-18.89 


2.19 


7.31 




9.54 


NGr 7fS83 


3.1 


5.2 


-19.53 


1.73 


7.63 




8.54 


NGr ?903 


4.0 


6.3 


-20.16 


1.38 


8.39 




9.01 


NGC ?Q7(S 


5.3 


1.6 


-17.35 


3.86 


6.42 


5 


7.49 


Nnr 303 i 


2.4 


2.4 


-19 90 


1 51 


7 4? 


u 


9 1 5 


NGr 3 1 84 


5.9 


7.7 


-19.11 


2.02 


8.35 




9.11 


NGr 3344 


4.0 


6.9 


-18.89 


2.19 


7.74 




9.01 


NGC 3351 


3.0 


8.3 


-19.46 


1.78 


8.08 




8.67 


NGC 3368 


1.8 


9.4 


-20.12 


1.40 


8.18 




8.95 


NGC 3486 


5.2 


8.2 


-18.84 


2.23 


7.50 




9.03 


NGC 3571 


4.0 


8.0 


-20.31 


1.30 


8.75 


5 


9.63 


NGr 3593 

. > V J V. JJ7J 


-0.4 


6.9 


-17.50 


3.65 


7.62 




7.75 


NGC 3623 


1.0 


8.9 


-20.17 


1.37 


7.62 




8.27 


NGC 3627 


3.0 


7.9 


-20.40 


1.26 


8.55 


6 


8.56 


NGC 3fS?8 


3.1 


9.2 


-20.67 


1.14 


8.62 




9.33 


NGr 40?o 


6.9 


9.2 


-17.31 


3.91 


6.60 


5 


8.05 


NGC 4062 

* iVIV T^V/VJZ. 


5.3 


9.4 


-18.78 


2.28 


7.63 




8.47 


NGr 409fS 


5.3 


7.9 


-19.49 


1.76 


7.75 


5 


8.86 


NGr 4 1 44 

. > V I V *T i. 1 1 


6.0 


3.1 


-15.93 


6.48 


6.31 




8.09 


NGr 4744 

l *i VIV tZ.t'T 


6.1 


2.3 


-18.06 


2.97 


6.62 




8.72 


NGr 4745 

. N VlV tZtJ 


0.1 


10.5 


-17.97 


3.07 


7.39 




6.61 


NGr 4274 


1.7 


10.9 


-19.33 


1.87 


8.27 




8.75 


NGr 4788 

i *i V I V ti,00 


7.1 


7.5 


-16.32 


5.62 


6.67 


5 


8.52 


NGr 43 1 4 

. > V I V tjlt 


1.0 


11.5 


-19.02 


2.09 


7.69 




6.43 


NGr 4359 


5.0 


14.3 


-17.49 


3.66 


6.55 




8.44 


NGr 4414 

! *i VIV 1 1 1 1 J. ""I" 


5.1 


8.9 


-19.25 


1.92 


8.48 




8.90 


NGr 4448 

! *i VIV 1 1 1 O 


1.8 


8.2 


-17.86 


3.20 


7.39 




7.38 


NGr 4490 

- iVIV I I /U 


7.0 


8.0 


-20.93 


1.04 


7.45 


6 


9.54 


NGr 4437 


6.0 


11.6 


-20.70 


1.13 


8.14 




7.90 


NGr 4575 


5.9 


13.5 


-18.11 


2.92 


6.57 




7.86 


NGr 4559 

* iVIV ^*J~>S 


6.0 


9.8 


-20.35 


1.28 


8.26 




9.57 


NGC 4565 

* i v I v i - ' vi - ' 


3.2 


13.8 


-21.74 


0.77 


8.62 




9.48 


NGr 4(S05 


4.9 


3.0 


-17.58 


3.54 


6.82 




8.05 


NGr 4(S3 1 


6.6 


7.9 


-21.46 


0.86 


8.03 




9.58 


NGr 47 3 (S 


2.4 


5.3 


-19.27 


1.91 


7.86 


5 


8.23 


NGr 48?(S 


2.4 


5.5 


-19.86 


1.54 


7.79 




8.07 


NGr 5055 


4.0 


7.3 


-20.43 


1.25 


8.80 


6 


9.40 


NGC 5 1 94 

. N VIV Jlyt 


4.0 


7.0 


-19.74 


1.61 


9.29 


5 


9.21 


NGr 5457 


5.9 


5.0 


-20.26 


1.33 


8.50 


5 


9.79 


NGr 6503 


5.9 


4.6 


-17.77 


3.31 


7.35 




8.86 


NGC 6946 


5.9 


4.1 


-20.12 


1.40 


8.74 




9.55 


NGr 7(S40 


5.3 


5.5 


-18.75 


2.31 


6.93 




9.62 


NGr 185 


-4.8 


0.7 


-13.83 


14.00 


4.81 


7 


5.18 


NGr 705 


-4.7 


0.7 


-13.61 


15.18 


4.95 


7 


5.57 




4 f. 


f, Q 

0.7 


1 fx 7^ 


5 81 


5 11 


7 


7 fD 


NGr 3265 


-4.8 


15.7 


-17.28 


3.96 


7.13 


7 


7.95 


NGC 3928 


-4.5 


12.1 


-17.35 


3.86 


7.36 


7 


8.22 


NGC 5128 


-2.1 


5.3 


-20.59 


1.17 


8.16 


7 


8.28 


NGC 5666 


6.4 


23.6 


-18.90 


2.19 


8.00 


7 


8.63 


NGC 1819 


-1.9 


44.8 


-20.23 


1.34 


9.10 


7 


9.13 


NGC 3032 


-1.8 


16.7 


-18.14 


2.89 


7.72 


7 


7.76 


NGC 4138 


-0.9 


10.9 


-17.97 


3.07 


7.13 


7 


8.54 


NGC 7465 


-1.9 


20.6 


-18.57 


2.47 


8.11 


7 


9.20 


NGC 3413 


-1.8 


7.9 


-16.66 


4.96 


7.21 


8 


7.95 
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Object 


T 


A/[/r'M P c] 


M B - 5 log h 


X 


log(M H2 /[X/r 2 M Q ]) 


Ref. H 2 


log(M ffl /|7i : 


NGC 5866 


-1.2 


9.5 


-19.23 


1.94 


7.81 


8 


8.15 


NGC 47 1 


-0.8 


13.8 


-19.02 


2.09 


8.25 


8 


7.20 


NGC 4459 


-1.4 


13.3 


-19.37 


1.84 


8.30 


8 


6.70 


NGC 4526 

. iVIV - T»JZ.V» 


-1.9 


6.7 


-18.63 


2.41 


8.30 


8 


7.05 


NGr 693 

» > V J V U7J 


0.1 


15.5 


-18.08 


2.95 


7.53 


8 


8.74 


NGr 7685 


-1.1 


11.0 


-18.32 


2.70 


7.45 


8 


8.79 


NGr 2773 


1.0 


20.7 


-19.47 


1.77 
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Figure Bl. Data points (subsample of the data shown in AppendixfA) rep- 
resent the observed relation between disc stellar mass and cold gas 
mass M g - ds . The solid line shows the best power-law fit and the shaded en- 
velope its l-o" uncertainty. This power-law is given in Eq. i23\ and has a 
slope of ff, = 1.46 ± 0.10. 

APPENDIX B: DIVERSE PHENOMENOLOGICAL 
RELATIONS 

This section summarizes the two phenomenological relations given 
in Eqs. l|23) and 

Bl Stellar mass versus gas mass 

From the galaxy sample presented in Appendix lAl we extracted all 
25 Scd/Sd-type galaxies (6 ^ T ^ 9), that is all objects approxi- 
mating pure discs. For these objects the total gas masses M gas were 
calculated via M gas = (Mm + M^)/{3. Additionally, we estimated 
the stellar mass M sta rs of each galaxy from the I-band magnitude M\ 
via dMo etal.lll998l) . 

log(M stilrs /M ) = 1.66 + log(Y,) - Mi/2.5, (Bl) 

where the mass/light - ratio log(Ti) = 1.2 has been adopted from 
iMcGaugh & de Blokl ( fl997h . 

The resulting data points displayed in Fig. IB II reveal an ap- 
proximate power-law relation between M gas and M stars . We have 
fitted the corresponding free parameters , y t to the data points by 
minimizing the x-y-weighted rms-deviations. l-<x errors for these 
parameters were obtained via a bootstrapping method that uses 10 3 
random half-sized subsamples of the 25 galaxies and determines 
the power-law parameters for every one of them. The standard 
deviations of the distributions for a\ and y\ are then divided by 
V2 to estimate 1-cr confidence intervals for the full data set. The 
best power-law fit and its 1-cr confidence interval are displayed in 
Fig. lB2l while explicit numerical values are given in Section [5~3l 

To first order, one would expect that M stars depends linearly on 
M gas , if both masses scale linearly with the mass of the parent haloe. 
The over-proportional growth of M slals (ai = 1.46±0.10> l)could 
be explained by the fact that more massive galaxies are generally 
older and therefore could convert a larger fraction of hydrogen gas 
into stars. 



Understanding the H2/HI Ratio in Galaxies 19 




Figure B2. Relation between disc scale radius r t ij sc and disc stellar mass 
AfjJjJ;. Squares represent 34 nearby spiral galaxies observed by Kregel et al. 
12002). Filled squares correspond to Sa/Sb-type galaxies, and empty 
squares represent Sc/Sd-type galaxies. These data were used to fit the free 
parameters of the model in Eq. 1241 . The solid line shows this model for 
T = 2 (i.e. Sab-type galaxies), while the dashed line shows the model 
for T = 6 (i.e. Scd-type galaxies). The slop e of these power-laws is 
Q'2 = 0.45 ± 0.05, consistent with a Freeman law (McGaugh et alj 1995t) . 



B2 Scale radius versus stellar mass 

lKregeletalJd2002h investigated a sample of 3 4 nearby edge-on spi- 
ral galaxies, drawn from the ESO-LV catalog l lLauberts & Valentiinl 
1989) using four selection criteria: (i) inclination i ^ 87 deg, (ii) 
blue diameter D^ s > 2.2 arcmin, (iii) Hubble type from SO-Sd, (iv) 
only regular field galaxies, i.e. no interacting systems, no warped 
or lopsided syste ms. This sample is complete in terms of sam- 
ple selection (see iKregel et alj|2002t lDavieslll990h . but the sam- 
ple volume i s too small to conta in rare objects. For each galaxy 
in the sample IKregel et al] f2002) determined the scale radius r disc 
of the stellar disc from the I-band luminosity profiles. They also 
obtained the morphological Hubble type T for each source from 
the Lyon/Meudon Extragalactic Database (LEDA). Additionally, 
we estimated the disc stellar masses Af$£ from the I-band mag- 
nitudes of the disc components according to Eq. (Bl\ . 

Using these data, we investigated the relations between MjJ c g , 
r disc and T. The data points shown in Fig. |B2] suggest the approxi- 
mate power-law with Hubble type correction of Eq. i24l . The best 
fitting parameters Qf2, 72, and 6 were obtained as in Section lBTI and 
explicit numerical values with errors are given in Section [53l 

To first order, the raise— relation can be understood in 
terms of a dark mat ter halo with an isothermal, singular, and spher- 
ical structure (e.g. iMo et ai] |l998h . This model predicts that the 
virial radius r v j, is proportional to the cubic root of the dark matter 
mass M DM at any fixed cosmic time. If r disc were proportional to 
r v j r and M* ar c s were proportional to Mqm, one would expect rdi SC to 
scale as (M^° s ) 1/3 . Our empirical result, cri = 0.45 ± 0.05, shows 
a slightly stronger scaling, consistent with the empirical Freeman 
law (q-2 = 0.5), according to whi ch disk galaxies have approxi- 
mately constant surface brightness dMcGaugh et alll995l) . 

The secondary dependence of the rdi sc -M* ar c s relation on the 
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Hubble type T probably has multiple reasons: (i) early-type galax- 
ies have more massive stellar bulges, which present an additional 
central potential that contracts the disc; (ii) bulges often form from 
disc instabilities, occurring preferably in systems with relatively 
low angular momentum, and hence early-type galaxies are biased 
towards smaller angular momenta and smaller scale radii; (iii) 
larger bulges, such as the ones of lenticular and elliptical galaxies, 
often arise from galaxy mergers, which t end to reduce the specifi c 
angular momenta and scale radii (see also lObreschkow et alj2009h . 



